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GUIDED RADIATION BEAMS IN FREE ELECTRON LASERS 


I. Introduction 


In many short wavelength free electron laser devices the radiation 
beam will not be confined or guided by a structure such as a waveguide. 
Furthermore, in order to provide high gain and efficiency, it is usually 
necessary for the interaction length (length of wiggler field) to be long 
compared to the diffraction length (Rayleigh length) associated with the 
radiation beam. In the FEL the tendency of the radiation beam to diffract 
away over a distance of a few Rayleigh lengths can be overcome by a 

focusing phenomenon arising from the resonant coupling of the radiation and 

1 2 

wiggler fields with the electron beam. ’ This radiation focusing effect 

3 

plays a central role in the practical utilization of the FEL. This 
phenomenon was first analyzed for the low gain FEL with transverse effects 
where it was shown that the diffractive spreading of the radiation beam 
could be overcome by a focusing effect arising from the modified index of 
refraction.* Experimental evidence indicating optical guiding in the FEL 
has also been observed recently. - ^ 

Optical guiding in FELs operating in the small signal exponential gain 

regime has been studied for the asymptotic behavior of the radiation 
8-11 

beam, indicating that it is possible for the propagation of self¬ 
similar transverse modes. Recently, a general formalism for optical 

focusing, guiding and steering, called the source dependent expansion (SDE) 

12 

Method, has been developed and applied to FELs. The SDE method is an 
excellent analytical and numerical technique for solving the wave equation 
that governs the FEL interaction. An envelope equation for the radiation 
beam in the FEL can be derived using the SDE method, and it is very 
appropriate for studying perfect guiding of the radiation beam in FELs 
operating in the exponential regime. We have obtained analytic expressions 
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for the spot size, wavefront curvature, phase shift and growth rate of the 
perfectly guided radiation beam in FELs operating in the Compton regime for 
different transverse profiles of the electron beam. The intrinsic 
efficiency of the FEL in the exponential gain regime with guided radiation 
can be calculated from these expressions, and from which the trapping 
potential and desired beam quality at injection can be estimated. It is 
found that high current rf linacs,^’*^ with their higher energy and better 
beam quality, are quite suitable for driving relatively short wavelength 
FELs to beyond saturation where the wiggler is tapered to enhance the 
efficiency. These results have been verified by simulations based on the 
SDE method for FELs driven by either induction or rf linacs. 

One of the consequences of optical guiding in FEL is the bending of 
the optical beam by a curved or misaligned electron beam.^ The SDE 
formalism allows us to obtain a condition on the curvature of the electron 
beam in an FEL that the radiation beam will remain guided. 












II. Refractive Index Associated with FELs 


In our model, the vector potential of an axially symmetric, linearly 


polarized, radiation field is 


A^(r,z,t) = A(r,z) e^/c-at).^ + c.c., 


where A(r,z) is the complex radiation field amplitude, w is the frequency 


and c.c. denotes the complex conjugate. 


The wave equation governing A D is 


r l?( r y * 5 2 /3z 2 -c- 2 3 2 /3t 2 A R . - ^ (2, 


where J (r,z,t) is the driving current density. Substituting Eq. (1) into 


Eq. (2) leads to the following reduced wave equation, 


[r 3r ( r 3r) + 2l c 3z a ^ r » 2 ^ = S(r,z,a), 


where a(r,z) = |e|A/m n c = |a|exp(i<{>) is the normalized complex radiation 


amplitude and we have assumed that a(r,z) is a slowly varying function of 


z, i.e., |(3a/3z)/a| « w/c. The source function, S, is given by, 


s *-rJ 


It is possible to relate the source function, S, to the index of 


refraction associated with the medium by noting that the wave equation 


for A R in a nonmagnetic, time-independent, nonlinear medium is 


_2 2 2 2 2 

(v - (n (r,z,a)/c )3 /3t )A R = 0, where n is the index of refractic 


associated with the medium and is, in general, complex and a nonlinear 


function of a(r,z). Comparing the reduced wave equation written in terms 


of n(r,z,a) with Eq. (3) we find that the source function can be written in 


terms of n, 


L_ 










( 5 ) 


S(r,z,a) = (»/c) 2 [l-n 2 (r,z,a)ja(r,z). 

The refractive index associated with FELs can be obtained from the 
following derivation, where a number of simplifying assumptions are made. 

For example: the beam electrons are monoenergetic without betatron 

15 

oscillations and that the radiation is of a single frequency. We write 
the nonlinear driving current density, J^ f as 

J x = -I e l n b ( r >v w (z)v oz js(z-z(t,t 0 ))dt 0 , (6) 

where n^(r) is the ambient beam density, Vq z is the axial electron velocity 
at z = 0 , tp is the time a given electron crosses the z = 0 plane, 
v w (z) = |e|A w /rm 0 c(exp(ik w z) + c.c.) e^/2, is the wiggle velocity, y is the 
Lorentz factor, A w is the vector potential amplitude of the planar wiggler 
field and k y =2n/X w is the wiggler wave number. Substituting Eq. ( 6 ) into 
the expression for S, Eq. (4) gives 


S 



5(t-x(z 0 ,t 0 ))/r, 


(7) 


where a y = |e|A w /mQC 2 , T = tg + J dz'/v z (z',tg) and the t^ integration is 

0 

over all entry times. Equating Eq. (7) with Eq. (5) and carrying out the 
integration over t^, we find the index of refraction associated with the 
FEL to be given by 

n fel (r,z,a) = 1 + (c^(r)/ 2 W 2 ) (7 ^ , ( 8 ) 

*0 


z 

where ♦ = J ( w/c + k y - iln(a/|a|) - w/v z (z,\|;q) jdz + <(<q is the relative 
0 




phase between the electron and the ponderomotive wave, = - cotg is the 
initial phase of a given electron and 

2n 

/V = (2JI} 1 J d* 0 is an ensemble average over the initial phases. The 
0 0 

radial profile of the index of refraction as given by Eq. (8) supports 
self-focusing of the radiation in an FEL. It should be noted, for 
completeness, that the relative phase satisfies the pendulum equation given 
by 

a 2 */dz 2 = 3k /3z - y" 2 («/c) f4'' 1 aa 2 /az -k a a sinJ. (9) 
v L w w w J 
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III. Radiation Beam Envelope Equation 


In order to solve Eq. (3) we will use the source dependent expansion 
12 

(SDE) method. This formalism has the merit that with only a few modes it 
permits an accurate solution of the wave equation throughout the 
interaction region. In this method, we choose the following representation 
for a(r,z) in terms of Laguerre-Gaussian functions, 


a(r,z) = Y, a m < z > L 


2 r 2 1 -<l-ia(z))r 2 /r 2 (z) 


~ ^ m' ' m 2. . ’ 

m tr (z)J 

where m = 0,1,2,•••. In Eq. (10), a m ( z ) are the complex amplitude 

coefficients, r g (z) is the radiation spot size, a(z) is related to the 

2 

radius of curvature of the radiation beam wavefront, R = - (»/2c)r s /a and 
L m is the Laguerre polynomial. Solving for the unknown quantities a m , r g 
and a in terms of the source term S allows us to completely describe the 
radiation dynamics. The representation in Eq. (10) is underspecified, 
since, when Eq. (10) is substituted into Eq. (3) and moments of the source 
function taken, there remain more unknown quantities than available 
equations. The additional degrees of freedom in our representation allow 
us to specify a particular functional relationship for the unknown 
quantities r g and a in such a way that, if the radiation beam profile 
remains approximately Gaussian, the number of modes needed to accurately 
describe the radiation beam is small. This yields the following first 
order coupled differential equations for r g and a, 


r' - 2ca/wr = - r H T , 
s s s I’ 















and a set of first order ordinary differential equations for the complex 
amplitudes a m (z), 

a' + A a m = - i |F - mBa m . - (m+l)B*a .1, (11c) 

m mm L m m -1 m+lj ' 

where H = F^/Sq, ' a 3/9z, and ( )^ j denotes the real and imaginary part 

of the enclosed function. In Eqs. (11), the functions A , B, and F are 

m’ ’ m 

given by 

A (z) = r'/r + i(2m+l) f(l + a 2 )c/to 2 - ar'/r + a'/2l, 
m ss v sss ) 

B(z) = - ^arVr s +(l-a 2 )c/wr 2 - a'/2j-i [rVr g - 2ac/wr 2 j, 

CD 

F m (z) = §5 J dC S(C,z)L m (Oexp(-(l + i a )C/ 2 ), 

0 

where C = 2 r 2 /r 2 . 

s 

Equations (11a) and (lib) can be combined to give the following 
envelope equation for the radiation beam 

r g + K 2 r s = 0 , ( 12 ) 

where 

K 2 =(2c /«) 2 ^-l+C 2 <sin^> 2 +2C<cos«J/>+(w/2c)r 2 C <sin\|/>jr~\ (13) 

C(z) = (2v/r)G(z)a w /|ap(z)|, measures the coupling between the radiation 

2 3 

and electron beam, v = («^q r^/ 20 ) = 1^/17x10 is Budker's constant, 1^ is 

2 2 

the electron beam current in amperes, G(z)=(l-f)/(1+f) and f(z)=(r b /r g ) 
is the filling factor associated with a Gaussian electron beam density 
profile. The first term on the right-hand side of Eq. (13) is the usual 
diffraction term, the second and third terms are focusing while the last 
term provides a focusing or defocusing contribution. In the high gain 






















IV. Guided Radiation Beams in the Exponential Gain Regime 


By considering the lowest order transverse mode (Gaussian profile) of 


the radiation beam, we find that the source term appropriate for the high 


gain Compton regime is, 


S(r,z) 


<Vr>/c> 2 (a w k v f B ) 2 


y(l+a 2 /2)(Ak-ir) 2 


a(r,z) ; 


(14a,b) 


where Ak and T are the wave number shift and growth rate respectively and 


fg is the usual difference of Bessel functions due to the linear wiggler. 


The lowest order mode is taken to have the form 


a(r,z) = aQ(0)exp(iJ(Ak-ir)dz'-(l-ia)r 2 /r 2 ). 


For the purposes of illustration, we will consider the Compton FEL regime 


in which the electron beam has a Gaussian density profile, ng(r) = n^exp 


(- r /r^). The conditions for a guided radiation beam require that the 


waist and curvature of the radiation beam remain constant, (r g = a' = 0). 


Setting r g = a'=0 in Eqs. (11a,b) and solving for T, Ak, r g , and a, (see 


Appendix), the following results for a guided beam are obtained. 


r = (l+a 2 ) - 1 (l+2f) - 1 r n , 


Ak=ocr, 


(16a,b) 


2,0x374 


X (1+a /2) ,1/4... 

w w ' f (1+ 


o 7 /4 , 1/2 

2 nrf B 


(l+3f/2)' 


« = (f/(2+3f))' 













where r o =2f B (v/y) 1/2 a w k w (l-*-a 2 /2)~ 1//2 and f=rjvr 2 is the filling factor. 


In the special case of f=l, 


r s (f=l) = 0.25 X wvy 


2 /0 .3/A 

'i'll/4 (1+a w /2) 

, 1/2 1/2 ‘ 
Yf B a w 


Similar procedures can be performed for other transverse profiles of 
the electron beam. Conditions for guided radiation beams in the Compton 
regime for Gaussian, parabolic, and flat-top transverse electron beam 
profiles are summarized in Table I. 

In Fig. 1, we show the spatial evolution of the radiation waist for 
the induction linac-driven FEL parameters in Table II. The parameters in 
Table II are consistent with Eqs. (16) and have been chosen to produce a 
guided radiation beam in the Compton exponential gain regime. The guided 
beam conditions can be shown to be stable,^ and Fig. 2 shows that 
irrespective of the initial value, the spot size asymptotes to the matched 
(guided) beam value. Figure 3 shows the evolution of the spot size for the 
rf linac-driven FEL parameters in Table III. As in Table II, the 
parameters in Table IV have been chosen to produce a guided radiation beam 
in the Compton exponential gain regime and are consistent with Eqs. (16). 

Free electron lasers driven by either induction or rf linacs could 
initially operate in the guided, exponential gain regime until saturation 
occurs. Immediately prior to saturation, the ponderomotive potential can 
be large enough, as in the above illustrations, to trap a significant 
fraction of the beam electrons. At this point, the wiggler field can be 
spatially tapered to achieve a significant increase in the operating 
efficiency and a somewhat smaller input signal into the FEL amplifier. 











To determine the viability of tapering the wiggler, prior to 
saturation, the trapping potential associated with the ponderomotive wave 
is needed. For linearly polarized waves, the fractional trapping potential 


|e|4> 


^ = 2/2 


rm 0 c 


ll+a /2j 


1/2 


(17) 


The radiation amplitude at saturation can be obtained from the intrinsic 
efficiency of the FEL. Using arguments based on electron trapping in the 
ponderomotive wave, we find that the intrinsic efficiency in the 
exponential (maximum) gain regime is 


n = Ak/k,. ( 18 ) 

Using the induction linac parameters in Table II as an illustration, we 
find that the intrinsic efficiency is h = Ak/k v = 0.66%. From this, the 
fractional trapping potential at the end of the exponential gain regime is 
lel^trap/troQ 0 ma kl n g it possible to trap the electron beam while 

tapering the wiggler field. In addition, the initial fractional energy 
spread of the electron beam must be somewhat less than f). This places a 
limitation on the fractional energy spread of the electron beam, SE/E^ < 

0.66%. One contribution to the beam energy spread is the transverse 

2 

emittance, 6E/E^ = (l/2)(e n /rj j ) . Therefore, the normalized beam emittance 

1 n 

must satisfy, e < (2Ak/k ) r, = 0.034 cm-rad. 
n w b 

Similar estimates can be carried out for the rf linac parameters in 
Table III. Even though the intrinsic efficiency is only 0.25%, the 
fractional trapping potential of 2% prior to saturation is still large 
enough to trap the electron beam and the wiggler field can be tapered. 
However, the small intrinsic efficiency puts a more stringent requirement 
on the beam quality, e < 0.007 cm-rad, in an rf linac-driven FEL. 









Figure 4 shows the relative power for ten transverse modes used in a 
simulation of the guided radiation beam for parameters in Table II. The 
fundamental mode is at least three orders of magnitude larger than any of 


the higher modes, indicating the SDE method is an excellent numerical 
scheme and the analytic results obtained with only the fundamental mode are 
well justified. 







Bending and Guiding of Radiation Beams 


Using the SDE formalism, it is possible to discuss the bending of a 
radiation beam by a curved electron beam in an FEL. For small 
displacements of the electron beam centroid, a nonaxisymmetric modal 
expansion similar to Eq. (10) can be performed and the spatial evolution of 
the centroid of the radiation beam found. Figure 5 shows the centroids of 
the electron and radiation beams for an FEL in the trapped particle regime 
with parameters given in Table II. Steering of the radiation beam by the 
electron beam is clearly demonstrated in this figure. 

It is interesting to consider the conditions under which the 
radiation beam could be guided by a curved electron beam. We denote the 
radial position by r = Rq + x, where Rq is the radius of curvature of the 
electron beam and x is the radial displacement from the center of the 
curved electron beam, as shown in Fig. 6. The FEL refractive index 
(correct to order x/Rq) is 

n = n fel + x/R 0’ ( 19 > 

where n^^ is given by Eq. (8) In the exponential gain regime, a guided 

radiation beam in a curved FEL is possible if R A > R . where 

v 0 ~ min 

R . = r /|Re(l-n, .) I. (20) 

mm s 1 fel 1 


Substituting the expressions for T, Ak and r g , from Eqs. (16), into (20) 
yields 











(21a) 



For a numerical example of R m j n » consider the following parameters, 

y * 100, 1=2 kA, r^ = 0.3 cm, a w = 1.72, f = 1 and f g = 0.85 (Table II). 

For these parameters, the minimum turning radius required for a guided 

radiation beam is R . = 455 m. 

min 









VI. Conclusion 


The source dependent expansion (SDE) method provides an excellent 
analytical and numerical technique for studying optical focusing, guiding 
and steering in FELs. We find that guided radiation beams in the FEL can 
be achieved in the Compton exponential gain regime but cannot be maintained 
in the high gain trapped particle (tapered viggler) regime. Conditions for 
guided radiation beam with different transverse profiles of the electron 
beam have been derived in the Compton exponential gain regime of an FEL. 

Free electron lasers driven by either induction linacs, such as the 
ATA, or high current rf linacs can operate in the guided, exponential gain 
regime until saturation occurs. At this point, the wiggler field could be 
spatially tapered so as to operate the FEL in the trapped particle regime 
in order to further increase the operating efficiency. 

We also examined the possibility of bending or steering radiation 
beams in FELs. We found a condition which places a lower limit on the 
radius of curvation of the electron beam necessary for the radiation to be 
guided along a curved path. 
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Table I 


Guided Radiation Beam Conditions for 
Compton Exponential Gain Regime 


ELECTRON 

BEAH 

PROFILE 





















Table II 


Parameters Associated with an Induction Linac-Driven 


FEL in the Exponential Gain Regime 


I b = 2kA, (v = 0.118) 

E fe = 50 MeV, (r = 100) 

r, = 0.3 cm 

b -3 

e < 34x10 cm-rad 
n 

Wiggler Field 

Wavelength .X = 8 cm 

w 

Wiggler Strength B y = 2.3 kG (a y = 1.72) 
Radiation Beam 

Wavelength X = 10.6 urn 

Spot Size r g = 0.25 cm, (Z^ = 2 m) 

(guided beam) 

e-folding length L g = 1/T = 94 cm 


Electron Beam 
Current 
Energy 
Radius 
Emittance 



n - Ak/k = 0.66% 
w 

P =660 MW (a = 7xl0 -4 ) 
sat 

l e l*trap /YTn 0 c2 = 6 -°* 







Table III 


Parameters Associated with an RF Linac-Driven 
FEL in the Exponential Gain Regime 


Electron Beam 
Peak. Current 
Energy 
Radius 
Emittance 

Wiggler Field (planar) 
Wavelength 
Wiggler Strength 


I fa = 500 A 
E fa = 150 MeV 

r, = 1 mm 

b -3 

e < 7x10 cm-rad 
n ~ 

A w = 12 cm 

B v = 900 G (a w = 1) 


Radiation Beam 

Wavelength A = 1 ym 

Spot Size r g (0) = 1.1 mm (Z R = 3.8 m) 

(guided beam) 

e-folding length L g = 1/T = 196 cm 


Intrinsic Efficiency 
Saturated Power 
Trapping Potential 


X] = Ak/k = 0.25% 
w 

P =180 MW (a = 7.25xl0~ 5 ) 
sat v 

l e l*trap /T V 2 ' n 
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Derivation of the conditions for guided radiation beam in the Compton 
exponential gain regime is shown in the following. When the transverse 
radiation beam profile is represented by the fundamental mode, Eq. (15), 
one can combine Eqs. (11a,b,c) to give the three dimensional dispersion 
relation for the Compton exponential regime, 

F F 

(Ak - in + ^-(1 - i«) + (-2 + -i) = 0 , (Al) 

<or a 0 a 0 

s 

where Fq, Fj are the overlap integrals of the source function, Eq. (14a), 
with the zeroth and first Gaussian-Laguerre modes. It can be shown that 

for the Compton source term, FQ=F^(l+2f), where f is the filling factor, 

2 2 

rj ) /r g . The dispersion relation Eq. (Al) is then reduced to 

9 F 1 

(Ak - iD + —(1 - ia) + 2—(1 + f) = 0 . (A2) 

»r a 0 

s 

By setting r g = 0 and a' = 0 in Eq. (11a,b) for guided radiation beams, we 
have, 

2 F 1 “ r 5 

(1- ia) 1 = -<-• (A3) 

a 0 c 

Evaluating F^/ag and substituting in Eq. (A3) gives the relations between 
the growth rate T, phase shift Ak and wavefront curvature factor a, 

Eqs. (16a,b). Substituting Eqs. (A3) into Eq. (A2) provides a second set 
of relations between T, Ak, a and the guided radiation beam radius r g , 






Ak + ^-[« 2 (1 + £)-£]- 0 , 


(A4) 


T + ^(2f -1) = 0 . 


Eliminating Ak and r from Eqs. (16a,b) and (A4) gives the relations of r g 
and a vith the filling factor f, Eqs. (16c,d). These algebraic equations 
can be easily solved numerically for the guided radiation beam values of 
r g , a, Ak and I*. 
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Fig. 5 Electron and radiation beam centroids, x b and x L for a displaced 
electron beam, x b = x c (l-sech(k c z)) with x c = r b /4 and 
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